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NOMENCLATURE 

a  =  3  Nondimensional  rectangular  channel  free  surface  length  (see  Fig.  2).  Dimensional 
length  equals  3 L. 

b  =  10  Half  length  of  the  nondimensional  rectangular  channel's  perfectly  conducting, 
moving  wall  (see  Fig.  2).  Dimensional  length  equals  10L. 

h  Nondimensional  current  density  stream  function  or  induced  magnetic  field  (see 
Eqs.  4.1  and  4.2) 

(*)  A  superscript  *  on  a  variable  denotes  a  dimensional  quantity. 

/*  Total  dimensional  load  current  across  the  channel 

h0  Nondimensional  parameter  representing  the  total  dimensional  current  current  /* 
where  /*  =  2 hQ  Ua  (crM)^2  (Ax*)  . 


Nondimensional  current  density  vector  = 


oU„  Br. 


Dimensional  height  of  channel. 
Nondimensional  pressure  (scalar)  = 


cUoBlL 


V  Nondimensional  fluid  velocity  =  j—  . 

U  o 

u  x-component  of  nondimensional  velocity. 

jy  ,  jz  Nondimensional  electric  current  density  components  in  the  y-  and  z-directions, 
respectively. 

* 

<p 


Nondimensional  electric  potential  = 


U„B0L 


<t )0  Nondimensional  electric  potential  between  moving  and  stationary  wall  of 

rectangular  channel. 

V  Nondimensional  gradient  (vector)  =  L  V  . 

(x*,  y*,  z*)  Dimensional  Cartesian  coordinates. 


(x,  y,  z)  Nondimensional  Cartesian  coordinates  =  ,  —  ,  ^-j,  (see  Fig.  2). 

(x,  v,  z)  Unit  vectors  in  Cartesian  coordinates  (see  Fig.  2). 

Ax*  Dimensional  length  along  rectangular  channel  =  L(Ax). 

B0  Dimensional  external,  homogeneous,  magnetic  field  =  B„  z  . 

cr  Fluid  electric  conductivity, 

p  Fluid  density. 

/.if  Fluid  viscosity. 

ob2[ 

N  Interaction  number,  a  nondimensional  parameter  =  — -y—  =  ratio  of  pondermotive 
force  to  inertial  force. 


Hartmann  number,  a  nondimensional  parameter  =  LBa  fr  =  positive  square 

v  Hf 

root  of  the  ratio  of  the  pondermotive  force  to  the  viscous  force. 

Magnetic  Reynolds  number,  a  nondimensional  parameter  =  U0  a  jJ.L  =  rado  of 
induced  magnetic  field  to  external  magnetic  field. 

Dimensional  velocity  of  perfectly  conducting,  moving  wall  of  rectangular  channel. 
Integral  expression  defined  as  a  residual  (see  Eq.  18). 

Total  nondimensional  Joulean  dissipation  per  unit  nondimensional  length  in 
rectangular  channel  (see  Eq.  24),  (Pf  =  HfUa2  Pj) 

Total  nondimensional  viscous  dissipation  per  unit  nondimensional  length  in 
rectangular  channel  (see  Eq.  23),  ( P ^  =  PfU02P ^ ) 

A  subscript  id  on  a  variable  denotes  the  one-dimensional  Couette  problem.  All 
variables  without  the  subscript  id  denote  the  two-dimensional  Couette  problem 
(see  section  entitled  "Comparison  of  Numerical  Parameters  for  One-dimensional 
and  Two-dimensional  Flow"). 


ABSTRACT 


Fully  developed,  viscous  liquid-metal  velocity  profiles  and  induced 
magnetic  field  contours  were  studied  for  Hartmann  numbers  ofM=2  and  10 
and  for  different  load  currents  for  a  particular  rectangular  channel 
configuration  ( two-dimensional  Couette  flow).  The  rectangular  channel 
was  assumed  to  have  a  homogeneous  external  (axial)  magnetic  field  parallel 
to  the  moving,  perfectly  conducting  top  wall  and  the  stationary,  perfectly 
conducting  bottom  wall.  The  two  stationary  side  walls  were  also  perfect 
conductors.  The  small  gap  between  the  moving  wall  and  each  side  wall  was 
an  insulating,  free  surface.  The  method  of  weighted  residuals  was  used  to 
obtain  truncated  series  solutions  for  the  variables  of  interest.  The  heavy 
load  currents  across  the  channel  were  obtained  by  simulating  an  external 
potential  to  the  conducting  moving  wall.  The  load  currents  in  each  case 
were  opposed  by  the  induced  electric  field.  Since  there  is  no  pressure 
gradient,  the  flow  along  the  channel  is  driven  by  the  viscous  effects  of  the 
moving  wall  and  the  Lorentz  body  force  and  is  retarded  by  the  stationary 
walls.  In  the  case  where  no  load  current  is  applied  across  the  channel,  the 
current  circulates  in  the  channel.  The  circulation  is  driven  by  the  generator 
that  is  due  to  the  axial  variation  of  velocity  in  an  axial  magnetic  field. 

The  numerical  data  presented  show  that  the  radial  gap  and  the  free 
surface  region  represent  electrical  resistances  in  parallel  between  the 
perfectly  conducting  stationary  wall  and  the  perfectly  conducting  moving 
wall.  The  data  also  show  that  the  resistance  of  the  radial  gap  increases  as 
M2  while  that  of  the  free  surface  increases  by  M  or  M1  ^  Thus,  as  M 
increases,  the  division  of  current  shifts  to  the  free  surface  region  and  the 
current  density  in  the  radial  gap  decreases  as  Mr1 .  The  theoretical 
magnetohydrodynamic  model  presented  here  was  developed  to  provide  data 
to  help  in  the  design  of  liquid-metal  current  collectors. 

Data  were  computed  for  one-dimensional  Couette  flow  with  no 
pressure  gradient  in  an  external,  homogeneous  axial  magnetic  field.  One¬ 
dimensional  Couette  flow  has  no  end  effects,  and  thus  the  data  were 
compared  with  corresponding  data  for  the  two-dimensional  Couette  flow 
case  to  determine  end  effects. 
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INTRODUCTION 


A  critical  technology  area  in  the  development  of  homopolar  generators  and  motors1'6 
is  the  mechanism  for  transporting  current  between  the  rotating  and  stationary  components 
of  the  system.  This  sliding  electric  contact,  or  current  collector  device,  must  successfully 
operate  in  the  very  high  current  density,  magnetic  fields  and  speed  regimes  of  a  homopolar 
system.  In  addition,  requirements  of  low  power  losses  and  long  life  time  make  traditional 
monolithic  brushes,  such  as  metal  graphite,  unattractive.  However,  liquid  metals  can 
provide  all  the  proper  operational  characteristics  as  well  as  minimizing  viscous  and 
electrical  power  losses  and  maximizing  the  system  lifetime.  For  instance,  state  of  the  art, 
continuous  duty  liquid  metal  current  collectors  have  successfully  transported  electric  current 
densities  of  3  x  107  amp/m2  in  model  devices.  Although  these  results  are  highly 
encouraging,  implementation  in  practical  systems  will  require  a  considerable  effort. 
Therefore,  further  improvements  in  current  collector  design  will  be  guided  by  a  thorough 
understanding  of  the  fundamental  principles  of  liquid  metal  flows  transporting  current  in 
high-magnitude  external  magnetic  fields.  To  this  end,  we  have  defined  in  this  paper  a 
magnetohydrodynamic  theoretical  problem  in  an  attempt  to  extract  some  of  the  basic 
physics  involved. 

In  contrast  to  the  many  studies7'9  of  fixed-wall  magnetohydrodynamic  duct  flow 
driven  principally  by  pressure  gradients  along  the  duct,  previous  work  by  the  authors 111-1 6 
has  been  directed  toward  solving  simple  models  of  collectors  as  two-dimensional 
rectangular  ducts  with  moving  walls,  transport  current,  and  magnetic  fields  of  varying 
orientation  driving  the  flow.10'15 

In  these  studies,  there  was  a  velocity  discontinuity  at  the  intersection  of  the  moving 
and  fixed  walls  which,  although  physically  impossible,  was  mathematically  tractible  and 
provided  a  number  of  fundamental  insights. 

The  purpose  of  the  present  paper  is  to  examine  a  more  physically  realistic  model, 
which  includes  stress-free  insulating  gaps  between  the  moving  and  fixed  walls.  As  with 
the  previous  models,  the  moving  and  fixed  conducting  walls  represent  the  rotor  and  stator 
of  a  current  collector  in  a  homopolar  device  w'hiie  the  two  stress-free  surfaces  represent  the 
spaces  in  between  the  rotor  and  stator  filled  with  a  liquid  metal  which  allows  device 
operation  with  relatively  low  friction  (see  Fig.  2).  The  problem  with  the  moving  perfectly 
conducting  wall  with  no  insulating  free  surfaces  in  the  comers  was  worked  by  the  authors 
and  served  as  a  basis  for  some  of  the  theoretical  work  herein  and  is  thus  presented  in 
Appendix  A.  The  magnetic  field  is  oriented  parallel  to  the  moving  wall  face,  a  geometry 
often  found  to  minimize  eddy  currents  in  practical  homopolar  generators.  The  parameters 
of  variation  are  applied  external  magnetic  field  strength,  fluid  electrical  conductivity, 
density,  viscosity,  gap  width,  rotor  axial  length,  axial  gap  width,  and  total  current  between 
stator  and  rotor.  Two-dimensional  results  are  presented  for  a  number  of  cases  for  channel 
fluid  velocity  profile,  induced  magnetic  field  in  the  channel,  and  Joulean  and  viscous  power 
losses,  and  these  results  an.  compared  with  similar  one-dimensional  solutions.  The  general 
analytical  solutions  for  the  fluid  velocity  and  current  density  can  serve  as  benchmarks  for 
various  numerical  computations  of  similar  problems.  The  values  of  parameters  chosen  lor 
numerical  calculations  are  typical  of  electrical  contact  parameters  in  generators  or  motors. 
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STATOR  LIQUID  METAL 


NOTE:  TOP  HALF  OF  COLLECTOR  SHOWN 


Fig.  1.  Typical  current  collector  configuration. 
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Fig.  2.  Rectangular  channel  cross  section. 
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NONDIMENSIONAL  MAGNETOHYDRODYNAMICS  EQUATIONS 


The  nondimensional  magnetohydrodynamics  equations  for  fully  developed, 
incompressible  laminar  duct  flow  in  an  external  axial  magnetic  field  can  be  expressed12'1? 
as 


A-1(v-v)v  =  -V/>  +  7  x?  +M'2V2V  , 

(1.1) 

^ 

/  =  -  Vtp  +  V  x  z  , 

(1.2) 

— * 

V  v  =  o  , 

(1.3) 

v-/  =  o . 

(1.4) 

gB20L 

In  Eq.  (1.1),  N  = -  represents  the  interaction  parameter,  the  ratio  of  pondennotive 

pUo 


force  to  inertial  force.  M  -  LB0 


is  the  Hartmann  number,  the  positive  square  root  of 


the  ratio  of  the  pondermotive  force  to  the  viscous  force.  It  is  assumed  during  this  work  that 
the  magnetic  Reynolds  number  Rm  =  U„  Of.iL  (ratio  of  induced  magnetic  field  to  external 
magnetic  field)  is  much  less  than  one.  The  nondimensional  variables  in  Eqs.  1.1  to  1.4  are 
defined  as 


V  =  fluid  velocity  normalized  by  U0 

— ►  * 

u0 

P  =  pressure  normalized  by  oU0BlL 
P* 

oU0BlL 

— ► 

j  =  electric  current  density  vector  normalized  by  oU0B0 

-** 

_  j 

oU0B0 

0  =  electric  potential  normalized  by  U0B0L 

i 

U0B0L 

V  =  gradient  normalized  by  L 
=  LV* 


(2.1) 

(2.2) 

(2.3) 

(2.4) 

(2.3) 


x,  y,  z  =  unit  vectors  in  Cartesian  coordinates  (see  Fig.  2).  (2.6) 


The  superscript  (*)  denotes  dimensional  variables.  The  x,  y,  and  z  coordinates  are 
normalized  by  L.  The  magnitude  B0  is  the  externally  applied  homogeneous  magnetic  field. 
U0  is  the  velocity  of  a  moving  wall  in  the  rectangular  channel.  The  variable  L  is  the  height 
of  the  channel.  The  characteristic  fluid  velocity  U0  is  the  velocity  of  the  moving  wall  in  the 
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rectangular  channel.  The  characteristic  length  L  is  the  height  of  the  channel.  The  liquid 
metal's  material  constants  are  electrical  conductivity  <7,  density  p  ,  magnetic  permeability  p 
and  viscosity  pp. 


RECTANGULAR  CHANNEL  PROBLEM 
MATHEMATICAL  SOLUTION 

The  problem  discussed  in  this  paper  consists  of  a  rectangular  channel  (see  Fig.  3) 
filled  with  a  liquid  metal  in  a  uniform  magnetic  field  B0  in  the  z-direction.  The  perfectly 
conducting  top  wall  is  assumed  to  move  with  velocity  u  =  U„.  The  moving  wall  extends 
from  z  =  -b  toz  =  +  /?atv  =  1.  At  each  top  comer  there  is  an  insulating  free  surface  in 
the  left  corner  from  z  =  -  (a  +  b)  to  -b  and  in  the  right  corner  from  z  =  +  h  to 
z  =  (a  +  b)  at  y  =  1.  The  perfectly  conducting  bottom  wall  is  stationary  with  velocity 
u  =  Q.  The  side  walls  are  also  perfect  conductors  with  velocity  u  =  0. 

The  flow  is  assumed  steady,  incompressible,  fully  developed  laminar  flow  in  the 

x-direction.  The  pressure  gradient  VP  is  set  equal  to  zero.  Thus,  there  is  no  pressure 

gradient  in  the  Jt-direction  along  the  channel  ^/dx~  0.  The  magnetic  field  is  uniform  and 

axial,  B0  =  B0  z.  No  assumptions  are  made  about  the  Hartmann  number  M.  The 
interaction  parameter  A  is  assumed  much  greater  than  one  in  this  problem.  By  symmetry 
considerations  numerical  results  in  the  paper  will  be  presented  only  for  the  right  half 
channel  0  <  z  <  (a  +  b)  with  0  <  y  <  1 . 

The  generalized  nondimensional  magnetohydrodynamics  equations  Eqs.  1.1  to  1.4 
can  be  transformed  into  the  following  vector  component  form 


M 


-2 


(~,2  n2 

o  u  d  u 

{dy2  dz2 

Jy  = - u  , 

3y 

,  _  30 

Jz  ’ 

dz 

djy  t  3 j2  _  q 

3  y  dz 


=  ~J> 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


Equation  3.1  is  the  ^-component  of  the  Navier  Stokes  equation.  Equations  3.2  and  3.3  are 
the  y-  and  z-components  for  the  current  density  Eq.  1.2  and  are  referred  to  as  Ohm's  law. 
Equation  3.4  comes  from  Eq.  1.4  and  is  referred  to  as  the  conservation  of  charge.  The 
continuity  equation  is  automatically  satisfied  in  Eq.  3.1  because  the  flow  is  uncompressible 

V  u  =  o) 
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Fig.  3  Rectangular  channel  with  moving  wall. 


These  equations  can  be  cast  in  terms  of  the  electric  current  density  stream  function  h, 
defined  by 


i  Bh 


(4.1) 


Jz  ~~ 


j  Bh 

M  By  ' 


(4.2) 


Replacing  the  current  densities  jy  and  jz  in  the  dimensionless  governing  equations  and 
eliminating  the  electric  potential  0  from  Ohm's  law  yields  the  following  set  of  second- 
order,  homogeneous  partial  differential  equations 


M 


M 


-,2 
B  U 

By2 

d2h 

By2 


B2u 


+  ■ 


Bz  2 
B2h 


Bz2 


Bz 


(5.1) 

(5.2) 


The  following  boundary  conditions  must  be  satisfied  by  these  equations  (see  Fig.  3). 
Due  to  the  symmetry  of  the  problem  only  the  right  half  of  the  rectangular  channel  is 
discussed.  The  same  boundary  conditions  apply  to  the  left  side  of  the  channel.  At  the 
moving  perfectly  conducting  wall  at  y  =  1  and  0  <  z  <  h  the  dimensionless  velocity  is  one 
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dh(l,z) 


the  free  surface  from  y  =  1  to  b  <  z  <  a  +  b,  the  shear  is  zero 


=  Oj.  Along 
and  the 


( u=  1)  and  the  z-component  of  the  current  density  jz  is  equal  to  zero 

\dy 

'MlA  =  0 

Wy  / 

current  density  stream  function  h(\,z)  equals  a  constant  h0.  Along  the  stationary  perfectly 
conducting  side  wall  from  0  <y  <  land  z  =  a  +  b  the  dimensionless  velocity  is  zero 

(dhlyya+b )  \ 

(u  =  0)  and  the  y-component  of  the  current  density  jy  equals  zero  - =  0  .  At  the 

\dz 

bottom  stationary  perfectly  conducting  wall  from  0 <z<a  +  b  at y  =  0  the 
dimensionless  velocity  u  is  zero  (u  =  0)  and  the  z-component  of  die  current  density  jz 


equals  zero  =  0 

lay 

at  z  =  0. 


du 

.  At  the  plane  of  symmetry,  —  =  0  and  h  =  0  from  0  <  y  <  1 

3z 


The  homogeneous,  second-order  partial  differential  Eqs.  5.1  and  5.2  are  assumed 
to  have  solutions  for  u(y,z)  and  h(y,z)  of  the  form 


'4y»z)=  X  U„(y)cos 

n~0 


n=0 


\l2n  +  n 

l  nz  \] 

[l  2  ) 

\a  +  bi\ 

\l2n+  11 

I JE.  |1 

LI  2  I 

(6.1) 


(6.2) 


which  satisfy  the  boundary  conditions  at  z  =  0  and  z  =  a  +  b. 

Substituting  series  expressions  6.1  and  6.2  for  u(y,z )  and  h(y,z)  respectively,  into  the 
simultaneous  partial  differential  Eqs.  5.1  and  5.2  and  solving  for  Hn(y)  results  in  the 
following  ordinary  differential  equation  for  Hn(y): 

-  2 Si  d^L  +  («f  +  S^]Hniy)  -  0  , 

dy 4  dy2 


(7) 


where 


\(2n  +  1 ' 

\  K 

2 

’( a  +  b )) 

Hn(y)  is  assumed  to  have  a  general  solution  of  the  form10 
H  (y)  =  coshPy  +  DsinhPy  . 

Substituting  Eq.  8  into  Eq.  7  gives  the  following  algebraic  equation  for  P 
PA  -  2SnP2  +  6n(Sn  +  M2)  =  0  , 
where  P\,i  =  V  S„  ±  iSnM  . 


(8) 

(9.1) 

(9.2) 
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Expressing  P\  and  P2  in  terms  of  real  and  imaginary  parts10,  we  have 


Pi  =  mr  +  imi 

(10.1) 

Pl-mr-  imi , 

(10.2) 

where 

/  'y  \  1  h 

mr  =  VI4[(l+M25n2)  +l]  , 

(10.3) 

and 

/  o\l h  ^  h 

mi  =  fZSn[{l+M28n)  -1 

(10.4) 

Thus  Hn(y)  has  the  following  general  solution 

=  >lni  cos^Piy)  +  cosh(P2y)  +  A„3  sinh(p\y)  +  4n4  sinhiPtf) , 

(11.1) 

or  more  simply 

Hn{y)  =  (A«i  +  An  2)  cosh(mry )  county) 

+  (Ani  +  An4)  sinh[mry)  cos(my) 

+i  {An3  -  An4)  cosh{mry)  sin(mjy) 

+i  (An  1  -  An 2)  sinh(mry)  sin(miy)  ,  (11.2) 

where 

@n\  ~  (Anl  +  An 2)  , 

@n2  ~  (An3  An 4)  , 

^«3  “  i  (An3  ~  An 4)  , 

an*  =  i(An\  ~Anl)  .  (11.3) 


Un(y)  can  be  determined  from  the  equation  (see  Eq.  5.2) 


1  \d2Hn(y) 

M  dy2 


SnHn{y) 


5nUn(y)=  0 


by  substituting  Eq.  11.1  for  Hniy).  Thus  Un(y)  is  determined  to  have  the  form 
Un(y)  =  An  1  icosh(P\y)  ~  An 2  icosh(P2y)  +  An3  isinf^Py)  -  An 4  isinl^P#') 
which  can  be  expressed  more  simply  as 

Un(y)  =  -  (-4„1  +  An 2)  sinh(mry)  sin(mjy) 

+  i  (An\  -  An2 )  cosh{mry )  co^my) 

+  i  (Ani  -  An4)  sinh[mry)  cos(my) 

-  i (An3  +  An4)  cosh(mry)  sin{my)  , 

where 

~(An  1  +  An  2)  =  -an  1  , 
i(^nl  ~An 2)  =  £Jn4  , 

i  {^n3  ~  ^n4)  =  Q/i3  < 

~{An3  An 4)  =  —  Qn 2  ■ 


02) 

(13.1) 

(13.2) 

(13.3) 
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When  the  boundary  conditions  at  y  =  0  are  applied ) 
Eqs.  6.1  and  6.2,  we  obtain 
mr 

&ri3  Qnl  and  <2^4  —  0  . 

mi 

Therefore,  the  expressions  for  Hn O')  and  Un(y)  are 
Hn{y)  =  an i  cos/t(mry)  cos(rmy) 

+  £/i2 


ldHn(y) 


\  dy 


=  0  and  Un{y)  =  0  to 


sinh(mry )  cos(miy)  -  —  cosh(mry)  sin(miy) 


and 


where 


Un(y)  =  -  an i  sinh(mry )  sin(miy) 

-  a,a\cosh{mry)  sin(mty )  +  ~  sinh(mry)  cos(mty 


mr  =  VI5fl[(l  +M2c5;2)1/2+  1 


wi  =  vr^[(i  +  m25„~2)1/2-  i 


1/2 

1/2 


and 


(14) 

(15.1) 

(15.2) 

(15.3) 

(15.4) 


(2n+\ \  K 
\  2  )f/i  + 


[a  +  b) 


Since  the  current  densities  jy  and  j2  are  even  and  odd  functions  of  z,  respectively,  h  is 
an  odd  function  of  z,  if  the  constant  value  of  A  at  2  =  0  is  chosen  as  zero.  The  boundary 
condition  jy-  0  at  the  free  surface  indicates  that  h(y,z)  is  a  constant  along  the  free  surface. 
The  constant  h0  denotes  this  value  and  represents  half  the  total  dimensionless  current  from 
the  stationary  walls  to  the  moving  wall. 

The  total  dimensionless  current  from  the  stationary  wall  to  the  moving  wall  along  the 
nondimensional  length  of  the  channel  Ax  is 
rb 


I  =  2(ax)  f  jy(\,z)dz  =  2(Ax)M  xh0  , 

Jo 


(16.1) 


where  the  dimensional  current  is  expressed  as 

/*  =  2h„  U0{opf)x  1/2  (Ax*)  ,  and  Ax*  =  t(  Ax*)  . 


(16.2) 


Once  we  obtain  the  computer  solutions  for  u(y,z)  and  h(y,z)  for  given  values  of  the  parameters 
a,  b,  M  and  h0,  we  can  compute  the  potential  difference  between  the  moving  wall  and  stationary 
walls  needed  to  drive  this  total  current. 


u(y,z )  +  A/-1 


dz 


dy 


for 


0  <z<b 


(17) 
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The  dimensional  potential  difference  is  then 

* 

<po  =  UoBoL  0O  . 

If  h0  has  the  same  value  for  a  range  of  values  of  B0,  it  represents  a  constant  total  electric 
current  from  the  stationary  wall  to  the  moving  wall  for  all  these  field  strengths.  On  the 
other  hand,  if  (J)0  has  the  same  value  for  a  range  of  values  of  B0,  it  represents  a  dimensional 
potential  difference  which  is  proportional  to  B0. 

A  positive  h0  represents  a  total  current  from  the  stationary,  perfectly  conducting 
bottom  wall  to  the  moving,  perfectly  conducting  top  wall.  A  negative  h0  represents  total 
current  flowing  in  the  opposite  direction.  Here,  we  assume  that  the  rotor  moves  in  the 
positive  x-direction  and  that  the  axial  magnetic  field  is  in  the  positive  2-direction.  In  a 
motor,  the  electromagnetic  (EM)  body  force  due  to  the  load  current  and  applied  magnetic 
field  is  in  the  same  direction  as  the  rotor's  motion.  Since  positive  h0  gives  an  EM  body 
force  in  the  positive  ^-direction,  results  for  positive  hQ  correspond  to  current  collectors  in 
motors.  In  a  generator,  the  EM  body  force  is  in  the  opposite  direction  to  the  rotor's 
motion,  which  is  the  case  for  negative  h0.  Therefore,  the  results  for  negative  h0 
correspond  to  current  collectors  in  generators. 


METHOD  OF  WEIGHTED  RESIDUALS 

From  the  boundary  conditions  at  the  moving  wall  and  at  the  free  surface,  we 
determine  the  coefficients  an\  and  an 2  by  defining  a  residual 

~  2  r 

sin(8nz)  +  1  -  X  U n(l) cos(8nz) 

L  n=0 


(18) 


where 


10 


and 


oo 


In  this  equation  we  have  used  the  notion  of  convergence  as  mean  square 
convergence.  The  first  integral  from  0  to  b  represents  the  convergence  of  the  appropriate 
series  to  the  boundary  conditions  along  the  moving,  perfectly  conducting  wall.  The  first 


1  ^/t  n 

squared  bracket  in  the  first  integral  represents  the  convergence  to  h  -  -  77  —  -  f . 

M  0y 

The  second  squared  bracket  in  the  second  integral  represents  the  convergence  of  the 
velocity  series  to  1  along  the  moving  wall  from  z  =  0  is  equal  to  h.  The  first  squared  term 
in  the  second  integral  represents  the  convergence  of  the  current  density  stream  function  or 
induced  magnetic  field  to  h0  along  the  free  surface  from  z  equals  b  to  a  +  b.  The  second 
squared  bracket  represents  the  convergence  of  the  v  partial  derivative  of  the  velocity  to  zero 
along  the  free  surface  from  b  to  a  +  b. 

The  residual  must  be  minimized  with  respect  to  the  coefficients 


dE 


=  0 


and 


da 


for  1  =  0,  1,2 


N  , 


(19) 


with  the  series  truncated  at  n  =/V.  The  result  is  2N  simultaneous  linear,  algebraic 
equations  with  the  2 N  unknowns 

ail,  «21.  -a,Vl  , 

012,  «22,  —dN2  ■  (20) 

For  a  specific  set  of  values  for  a,  b,  M ,  h0,  and  N,  we  solve  the  simultaneous 
equations  for  the  2 N  unknown  coefficients.  The  series  solutions  for  u(y,z)  and  h(y,z) 
now  give  the  values  of  u(y,z)  and  /i(y,z)  throughout  the  liquid  metal.  In  addition,  the 
following  relationship  holds  for  the  electric  potential 

<p{y,z)  =  -A/-1  £  dH-r'^  -  8nl  mv(5„z)  .  (21) 

nA)  - 


RESULTS  AND  DISCUSSION 

Graphical  data  are  presented  for  the  nondimensional  velocity  u(y,z)  and 
nondimensional  current  density  stream  function  or  induced  magnetic  field  h{y,z)  for  the 
magnetohydrodynamic  channel  with  a  moving  perfectly  conducting  top  wall  having  two 
free  surfaces  at  the  top  comers  (see  Fig.  4).  The  channel  parameters  were  chosen  to 
illustrate  variations  of  several  key  features  that  are  typical  for  current  collectors  in  motors 
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Fig.  4.  Coordinates  for  Contour  plots. 
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and  generators.  The  radial  gap  was  chosen  as  length  L  ,  each  axial  free  surface  gap  as  3L, 
and  the  length  of  the  moving  wall  as  20 L.  The  stationary  walls  and  moving  wMIs  were 
chosen  to  be  perfect  conductors.  These  lengths  were  nondimensionalized  with  respect  to 
L.  The  applied  magnetic  field  was  varied  through  the  Hartmann  number  M  =  B0L  *\J 

Typical  values  for  liquid  metal  current  collectors  in  homopolar  motors  would  be  M  =  2  to 
50.  Although  in  the  discussion  that  follows  the  dependence  of  various  quantities  on 
variations  of  the  Hartmann  number  has  been  interpreted  in  terms  of  applied  magnetic  field 
variations,  the  discussion  is  equally  applicable  to  constant  applied  magnetic  field  and 
variations  of  other  parameters  that  make  up  the  Hartmann  number.  These  include  the 
distance  between  the  top  and  bottom  conducting  walls  and  the  fluid's  electrical  conductivity 
and  viscosity.  In  this  work  the  nondimensional  pressure  gradient  was  set  equal  to  zero. 
Thus  the  fluid  motion  is  due  only  to  the  drag  of  the  moving  wall  and  the  j  xB0  body  force. 
Only  the  right  half  of  the  channel  was  drawn  for  each  figure  presented  herein  since  the 
calculated  value  of  the  velocity  is  symmetrical  and  since  the  induced  magnetic  field  is  anti- 
symmetrical  with  respect  to  z. 

The  total  current  from  the  stationary  to  the  moving  wall  (stator  to  rotor)  along  the 
length  of  the  channel  x  is  given  by 

I*  =  2{A  x*)oU0B0Lh0M-x  (22.1) 

I*  =  2/zc  Ua(oHf  )1^2  (a**)  ,  (22.2) 

where 

/  *  =  total  dimensional  current, 

U0  =  velocity  of  moving  wall, 

o  =  conductivity  of  fluid, 

fif  =  viscosity  of  fluid, 

hQ  =  parameter  representing  the  total  current, 

(Ax*)  =  (Ax)  L  =  dimensional  longitudinal  length  along  channel. 

The  total  nondimensional  viscous  dissipation  per  unit  nondimensional  length  in  the 
channel  is  calculated  from 

P,= 


where  Pjl  = 


r(a+b) 


-(a+b)  J  0 


8u 

U 


\2 


jdu 


dydz  , 


(23) 


HfVl  and  the  superscript  (*)  represents  a  dimensional  quantity.  Also,  the 
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total  nondimensional  Joulean  dissipation  per  nondimensional  unit  length  in  the  channel  is 
calculated  from 


!*( a+b )  r\ 


Pj  =  M 2 


[jy  +  jl]  dy  dz  , 


- (a+b)  J 0 


(24) 


where  P*  =  p/Ul  Pj  . 


Data  were  calculated  numerically  on  the  computer  for  half  the  nondimensional  length 
of  the  moving  wall  b  - 10  and  the  free  surface  nondimensional  length  a  =  3.  The 
Hartmann  number  cases  studied  were  M  =  2  and  10.  The  total  current  between  the  moving 
and  stationary  walls  was  represented  by  h0  =  0,  ±  20,  ±  40,  ±  60,  ±  80,  ±100.  A 
positive  h0  represents  the  total  current  moving  from  the  stationary  conducting  wall  to  the 
moving  conducting  wall,  and  a  negative  h0  represents  current  from  the  moving  conducting 
wall  to  the  stationary  conducting  wall.  The  quantity  h0  was  chosen  to  represent  typical  load 
currents  carried  by  current  collectors  in  generators  or  motors.  The  higher  positive  and 
negative  h0  values  correspond  to  very  high  load  currents.  Due  to  symmetry  considerations 
only  one  half  of  the  channel  is  shown  (i.e.,  right  side  of  channel).  The  y-axis  has  been 
stretched  by  a  factor  of  5  in  the  contour  plots  0  <y  <  1  and  the  z-axis  runs  from  0  to  13. 

Due  to  the  contour  routine  utilized  to  plot  the  contours,  the  singularity  at  the  top 
comer  of  the  moving  wall  has  been  smeared  across  two  elements,  one  of  which  lies  on  the 
free  surface.  Therefore,  current  lines  that  appear  to  enter  the  free  surface  actually  enter  the 
moving  wall,  so  that  current  densities  at  the  moving  wall  are  larger  than  those  shown  in  the 
contour  plots. 

Figure  5a  presents  composite  contour  plots  of  nondimensional  velocity  w(y,z),  and 
contour  plots  of  current  density  stream  function  (induced  magnetic  field)  h(y,z)  x  1000  in 
the  right  hand  side  of  the  channel  at  Hartmann  number  M  =  2  and  h0  =  0.  The  velocity 
contours  representing  flow  in  the  positive  x-direction  are  represented  by  solid  lines.  Later 
in  the  report,  velocity  contours  representing  flow  in  the  negative  x-direction  are  represented 
by  dashed  lines.  The  induced  magnetic  Field  contour  lines  are  represented  by  dashed  lines 
with  dots  and  arrows  showing  the  direction  of  current  density  flow.  When  h0  =  0  there  is 
no  net  current  flow  between  the  moving,  perfectly  conducting  top  wall  and  the  stationary, 
perfectly  conducting  bottom  wall.  The  analogous  one-dimensional  Couette  flow  problem 
(Appendix  B)  in  a  constant,  homogeneous  axial  magnetic  Field  with  no  nondimensional 
pressure  gradient  would  give  the  nondimensional  velocity  uiy)  =y,  the  external  applied 
nondimensional  potential  (p0  =  -0.5,  and  the  y-component  of  the  nondimensional  current 
density  jy  =  0.  In  Fig.  5a  for  the  region  0  <  z  <  9  and  0  <  y  <  1  we  have 
approximately  the  simple  Couette  flow  solution  w(y)  =y.  For  the  region  10<z<  13  and 
0<y  <  1  the  stress-free  surface  does  not  drive  the  fluid,  and  the  stationary,  perfectly 
conducting  sidewall  retards  the  fluid.  The  fluid  that  flows  in  this  region  moves  because  of 
its  viscous  interaction  with  the  simple  Couette  flow  between  the  moving  top  wall  and 
stationary  bottom  wall.  The  flow  pattern  extends  from  approximately  z  <  9  to  z  <  1 2. 
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The  induced  magnetic  field  acts  as  a  stream  function  for  the  current  density  and 
therefore  expresses  the  fact  that  the  direction  of  the  tangent  to  a  streamline  is  the  direction  of 
the  current  density.  The  current  density  at  each  point  of  the  fluid  is  determined  by  the 

Ohm's  law  relationship  j  =-V<J  +  iT  xz.  An  axial  variation  of  u(y,z)  acts  as  a  generator 

in  an  axial  magnetic  field.  This  will  drive  an  external  circulation  of  current.  In  the  figure 
there  is  a  large  axial  velocity  gradient  between  z  =  9.5  and  z  =  11.5.  Since  the  flow  is  in 
the  positive  x-direction,  the  current  density  circulation  is  counterclockwise.  As  the  current 
density  circulates,  the  Lorentz  body  force  tries  to  minimize  the  axial  variation  of  the 
nondimensional  velocity  u(y,z). 

The  axial  lines  of  the  current  density  circulation  are  inside  the  perfect  contours 
because  this  represents  a  lower  resistance  than  the  liquid  metal.  There  are  downward 
current  density  lines  for  0<z  <9.5  and  upward  current  density  lines  from 
9.5  <z<  13.  The  upward  current  density  lines  must  converge  to  the  edge  of  the  perfectly 
conducting  moving  wall  aty  =  1,  z  =  10  because  of  the  insulating  electrical  properties  of 
the  free  surface. 


If  the  problem  is  viewed  in  another  way,  the  externally  applied  potential 
<t>o  =  <Hy  =  1  ,z)  for  -b  <  z  <  b  must  produce  an  electric  field  at  each  point  of  the  fluid  that 
balances  some  global  average  of  if  x  f  throughout  the  fluid  so  that  there  is  no  net  current. 
If  there  were  pure  Couette  flow,  then  <p0  =  -0.5  (see  Appendix  B).  However,  the  fluid 
near  z  =  13  is  moving  at  u  <  y,  so  1 0O\ <  0.5.  For  0  <  z  <  9.5,  the  flow  approximates 


pure  Couette  flow;  thus  u  x  z  =~uy  is  larger  in  magnitude  than 


>  0,  so  current 


density  lines  are  downward.  For  z  >  9.5,  u  x  z  is  very  small,  so  Ey  >  u  and  current 
density  lines  are  upward.  The  current  density  decelerates  flow  for  0  <  z  <  9.5  and 
accelerates  flow  for  9.5  <z  <  13,  but  for  M  -  2  the  Lorentz  force  is  too  small  in 
magnitude  to  significantly  alter  u(y,z). 


Figure  5b  presents  composite  contour  plots  of  nondimensional  velocity  u(y,z)  and 
contour  plots  of  current  density  stream  function  (induced  magnetic  field)  h(y,z)  x  1000  in 
the  right  hand  side  of  the  channel  at  Hartmann  number  M  =  10  and  h0  -  0.  The  Lorentz 
body  force  acting  on  the  fluid  for  M  =  10,  h0  =  0  is  much  greater  than  for  M  =  2,  h0  =  0 
(Figs.  5a  and  5b).  For  the  case  where  M  =  10,  the  flow  profile  is  altered  and  the  axial 
velocity  gradient  is  reduced,  in  contrast  to  the  M  =  2  case.  Figure  5b  indicates  that  the 
nondimensional  velocity  profile  for  9.5  <z  <  13  has  been  increased  slightly  relative  to 
that  for  M  =  2.  We  need  to  remember  that  the  vertical  scale  is  stretched  by  a  factor  of  5. 
The  effect  of  the  moving  wall  on  the  fluid  flow  is  extended  well  beyond  the  end  of  the 
moving  wall.  Correspondingly  the  velocity  for  0  <  z  <  9  has  been  reduced  below  that  for 
Couette  flow.  For  the  current  density  stream  function  there  are  several  important  features 
to  recognize.  For  M  =  2,  ho  =  0  and  9.5  <  z  <  13,  virtually  all  of  the  current  density  lines 
move  from  the  stationary  wall  at  y=0  and  converge  to  the  end  of  the  moving  wall.  For 
M=  10,  a  large  number  of  current  density  lines  leave  the  sidewall  at  z  =  13  for 
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Fig.  5b.  M  =  10,  h0  =  0 

Fig.  5.  Nondimensional  velocity  contours  u(y,z)  and  nondimensional  current  density 
stream  function  h(y,z)  x  1000  for  Hartmann  numbers  M  =2  and  10,  and 
h0  =  0. 
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0.25  <z<  1.  There  are  several  ways  to  view  this.  The  solution  for  one-dimensional 
Couette  flow  (presented  in  Appendix  B)  showed  that  a  current  across  the  magnetic  field 
lines  produces  a  velocity  change  which  in  tum  opposes  the  current  density.  This  leads  to 
an  effective  resistance  to  the  current  density  jy  which  is  given  by 

R*ff=R[l  +!Jfj  •  (25) 

On  the  other  hand,  current  density  lines  along  magnetic  field  lines  produce  no  body  force 
j  x  B0,  so  these  lines  do  not  alter  the  flow  or  generate  any  increase  in  the  effective 

resistance.  Therefore,  jz  always  sees  a  resistance  R ,  while  jy  sees  the  effective  resistance 
R  =  Reff.  For  the  case  M  -  2,  h0  =  0,  there  is  a  small  difference  between  the  resistances 
for  jy  and  j2,  so  both  current  densities  see  approximately  the  same  resistance  in  the  fluid. 
The  distance  from  the  stationary,  perfectly  conducting  wall  at  y  =  0  to  the  moving, 
perfectly  conducting  wall  at  y  =  1  is  of  length  L,  but  the  approximate  length  from  the 
stationary,  perfectly  conducting  side  wall  at  z  =  13  to  the  moving  wall  is  3 L,  so  the  current 
density  path  of  least  resistance  is  from  the  bottom  to  the  moving  wall.  For  the  case 
M  =  10,  the  effective  resistance  Rejj  =  9.33 R,  so  that  the  path  from  the  sidewall  to  the 
moving  wall  has  a  lower  resistance  than  that  from  the  stationary  bottom  wall  to  the  mov  ing 
top  wall.  Of  course,  the  current  density  cannot  be  only  the  jz  -component,  because  as  it 
reaches  the  moving  wall  it  must  tum  into  the  jy  -component  to  enter  the  moving  wall. 
However,  the  current  density  lines  have  only  a  short  distance  to  travel  in  the  y-direction. 

An  M~x V-  thick  free  surface  boundary  layer  would  predict  that  the  current  lines  would 
be  a  series  of  parabolas  converging  to  a  current  sink  at  the  edge  of  the  moving  wall  and  the 
current  density  lines  in  Fig.  5b  have  precisely  this  parabolic  shape.  The  M~1/2  layer  would 
have  an  appropriate  actual  thickness  of  where  H  is  the  dimensionless  height  of 

the  layer  along  a  magnetic  line. 

A  second  feature  of  Fig.  5b  is  the  recirculation  of  the  current  density  lines  near 
z  =  10.5,  y  =  0.  In  spite  of  the  Lorentz  body  force  effects,  the  fluid  near  z  =  13,  y  =  0 
remains  at  small  velocity  because  of  the  no-slip  condition  on  the  two  walls  at  the  corner. 
Therefore,  there  is  a  local  battery  action  at  this  comer.  One  part  of  the  circulation  is  short- 
circuited  through  the  stationary  bottom  wall,  but  the  other  part  of  the  path  closes  through 
the  liquid.  There  would  be  a  higher  resistance  in  a  path  to  the  perfectly  conducting  moving 
wall  than  in  the  path  just  inside  the  liquid. 

Figures  6a  through  6c  present  composite  contour  plots  of  nondimensional  velocity 
u(y\z)  and  current  density  stream  function  or  induced  magnetic  field  h(y,z)  for  Hartmann 
number  M  =  2  and  h0  =  20,  60  and  100.  The  stream  function  contours  are  represented  by 
dotted  lines  with  arrows.  The  velocity  profile  contours  in  the  positive  .x-direction  are 
represented  by  solid  lines.  The  positive  values  of  h0  represent  relatively  strong  net  current 
moving  from  the  stationary,  perfectly  conducting  bottom  wall  to  the  moving,  perfectly 
conducting  top  wall.  In  each  h0  case  for  M  -  2  there  are  essentially  two  current  paths 
available.  In  the  first  path  current  can  flow  across  the  relatively  narrow  gap  between  the 
top  moving  wall  and  the  bottom  stationary  wall.  However,  the  large  positive  y-component 
of  the  nondimensional  current  density  jy  accelerates  the  flow,  and  then  the  induced  electric 
field  opposes  the  current,  increasing  the  effective  resistance  by  a  small  amount.  The 
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second  current  path  is  along  the  magnetic  field  lines  near  the  free  surface  and  then  into  the 
edge  of  the  moving  wall.  This  has  a  longer  path  than  the  first  process,  approximately  3 L. 
The  current  density  lines  along  the  field  lines  do  not  produce  a  flow  acceleration,  thus  this 
phenomenon  decreases  the  effective  resistance.  However,  the  final  parts  of  the  current 
density  paths  involve  a  very  large  nondimensional  current  density  jy  near  the  edge  of  the 
top,  moving  wall  which  locally  accelerates  the  flow  and  which  produces  an  increased 
effective  resistance  for  this  small  region.  The  effective  resistance  of  the  end  of  this  moving 
wall  end  region  probably  varies  as  some  power  of  M W  or  M.  The  resistance  across  the 
gap  increases  as  M2.  Therefore,  an  increase  of  M  will  cause  a  progressive  shift  of  the 
current  from  the  radial  gap  of  length  L  to  the  region  near  the  free  surface.  The  larger  the 
transport  current  the  faster  the  fluid  moves  in  the  channel  if  the  other  parameters  stay  the 
same. 

The  nondimensional  positive  velocity  profiles  u(y,z)  for  the  two-dimensional  model 
when  M  =  2  in  the  region  0  <  z  <  10  and  0  <y  <  1  approximate,  but  are  slightly  less 
than  the  corresponding  nondimensional,  one-dimensional  velocity  profiles  for  the  Couette 
model  for  M  =  2  and  h0  =  20,  60,  and  100.  Appropriate  values  of  the  nondimensional 
velocity  u(y)  -  0.25,  0.50,  and  0.75)  and  u^x  for  one-dimensional  Couette  flow  with  an 

axial  magnetic  field  and  no  pressure  gradient  dP/fix  are  presented  in  Table  1  for  M  =  2  for 

comparison.  The  nondimensional  velocity  profiles  are  similar  for  the  two-dimensional  and 
one-dimensional  model  because  current  end  effects  are  small  in  the  two-dimensional  model 
at  a  low  Hartmann  number  of  M  =  2.  That  is,  all  the  two-dimensional  figures  for  h0  =  20 
to  100  show  all  induced  current  density  contour  lines  in  the  region  0  <  z  <  10.5  for 
M  -  2  from  the  stationary,  perfectly  conducting  wall  to  the  moving,  perfectly  conducting 
wall  except  for  one  current  density  contour  across  the  free  surface  shown  in  all  figures 
except  Fig.  6c.  In  all  figures  for  M  =  2,  as  positive  h0  increases,  the  induced  magnetic 
field  contour  lines  generally  increase  in  magnitude  and  the  nondimensional  velocity  contour 
lines  also  generally  increase  in  magnitude. 

Figures  7a  through  7c  represent  the  composite  contour  plots  of  nondimensional 
velocity  u(y,z )  in  the  positive  x-direction  and  nondimensional  induced  magnetic  field  for 
M  =  10  and  h0  =  20,  60,  and  100.  The  positive  values  of  h0  represent  relatively  strong 
net  current  between  the  stationary,  perfectly  conducting  bottom  wall  and  the  moving, 
perfectly  conducting  top  wall.  The  two-dimensional,  nondimensional  velocity  profiles 
u(y,z)  for  the  region  0  <  z  <  10  for  0  <y  <  1,  are  in  general  significantly  less  than  the 
corresponding  one-dimensional,  nondimensional  velocity  profiles  u(y)  of  simple  Couette 
flow  in  an  axial  magnetic  field  for  M  -  10  and  h0  =  20  to  100.  Appropriate  values  of  the 
nondimensional  velocity  u(y=  .25,  .50,  and  .75)  and  umax  for  one-dimensional  Couette 
flow  for  M-  10  are  shown  in  Table  1.  As  the  Hartmann  number  M  was  increased  from  2 
to  10  in  this  set  of  computed  data,  the  axial  magnetic  field  was  increased,  which  results  in 
the  figures  for  M  =  10  showing  many  more  induced  magnetic  field  contour  lines  than  the 
figures  for  M=  2  from  the  perfectly  conducting,  stationary  right  side  wall  to  the  tip  of  the 
perfectly  conducting  moving  wall.  The  induced  magnetic  field  contour  lines  for  M  =  2  for 
the  region  0<z<  10  and  0<y<  1  are  the  same  magnitude  as  the  corresponding 
contour  lines  for  M=  10,  except  that  in  the  M  =  2  case  the  lines  are  sometimes  shifted 
slightly  to  the  left.  Thus,  there  is  a  slight  decrease  of  current  across  the  gap  of  length  L 
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Table  1. 

One-dimensional 

Couette  flow 

values  for 

axial  magnetic 

field. 

Hartmann  No. 

h0 

«(•  75) 

u(.  50) 

u(.  25) 

Umax 

HHMi 

20 

1.13 

1.00 

1.50 

0.63 

1.13 

0.75 

40 

1.51 

1.01 

1.57 

0.63 

M  =  2 

60 

1.89 

2.00 

1.39 

2.04 

0.58 

80 

2.27 

2.50 

1.77 

2.53 

0.56 

100 

2.65 

3.00 

2.15 

3.03 

0.55 

20 

2.65 

3.00 

2.15 

3.03 

0.55 

40 

4.55 

5.50 

4.05 

5.51 

0.53 

M=  10 

60 

6.40 

8.00 

5.95 

8.01 

0.52 

80 

8.35 

10.50 

7.85 

10.50 

0.51 

100 

10.25 

13.00 

9.75 

13.05 

0.51 
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from  the  bottom  conducting  wall  to  the  top  moving  conducting  wall  when  M  goes  from  2  to 
10  as  would  be  expected  from  previous  arguments.  The  nondimensional  velocity  profiles 
for  the  Af  =  10  case  are  generally  of  a  much  higher  order  of  magnitude  than  the 
corresponding  nondimensional  velocity  profiles  for  the  M  =  2  case  because  of  the  much 
increased  j  x  B0  body  force  in  the  fluid.  The  opposing  induced  electric  field  a(v  x  B0) 
is  much  larger  in  the  M  =  10  case  than  the  M  =  2  case. 

These  sets  of  data  for  M  -  2  and  M  =  10  for  hQ  =  20,  60,  and  100  support  the 
concept  that  the  radial  gap  and  the  free  surface  region  represent  electrical  resistances  in 
parallel  between  the  stationary  walls  and  the  moving  wall.  It  also  shows  that  the  resistance 
of  the  radial  gap  increases  as  A/2  while  that  of  the  free  surface  region  increases  as  M  or 
A/!/2.  Thus  as  M  increases,  the  division  of  current  shifts  to  the  free  surface  region  and  the 
current  density  in  the  radial  gap  decreases  as  M_1-  We  have  already  justified  the  idea  that 
the  radial  gap  resistance  varies  as  M  2.  The  resistance  of  the  free  surface  region  is  less 
clear.  In  the  limit  as  M co,  this  region  becomes  a  parabolic  one.  Figure  8  shows  high 
Hartmann  number  regions  at  tip  of  the  moving  wall.  All  current  must  pass  through  the 
O  (AT1)  x  O  (AT1)  comer  region  at  the  edge  of  the  moving  wall.  However,  from  this 
theoretical  development,  it  is  not  clear  whether  the  free  surface  layer  with  A/'1/2  thickness 
or  the  comer  region  with  A/'1  thickness  is  the  dominant  resistance. 


Fig.  8.  High  Hartmann  number  regions  at  tip  of  moving  wall. 


Figures  9a  through  9c  present  the  composite  contour  plots  of  nondimensional  velocity 
u(y,z)  and  nondimensional  induced  magnetic  field  h(y,z)  for  M  =  2  and  ha  -  -  20,  -  60, 
and  -  100.  The  dashed  lines  with  dots  show  velocity  profiles  in  the  negative  direction. 

The  solid  lines  are  velocity  profiles  in  the  positive  ^-direction  as  in  previous  contour  plots. 
The  induced  magnetic  field  contour  plots  are  dashed  lines  with  arrows.  The  negative 
values  of  h0  represent  relatively  strong  net  current  from  the  moving,  perfectly  conducting 
top  wall  to  the  stationary,  perfectly  conducting  bottom  wall.  In  each  negative  h0  case  there 
are  two  current  paths  available  as  discussed  previously.  In  the  first  case  current  can  flow 
directly  from  the  moving  wall  to  the  stationary  wall.  However,  the  large  negative  v~ 
component  of  the  nondimensional  current  density  jy  accelerates  the  flow  in  the  negative  v- 
direction  and  the  induced  electric  field  opposes  the  current,  increasing  the  effective 
resistance  across  the  channel  by  a  small  amount.  The  velocity  of  the  fluid  is  in  the  positive 
x-direction  close  to  the  moving  wall  because  of  viscous  action  between  the  moving  wall 
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and  the  fluid.  The  second  current  path  is  out  of  the  edge  of  the  moving  wall  and  along  the 
magnetic  field  lines  near  the  free  surface.  As  previously  discussed,  a  shift  of  M  from 
M  =  2  to  M=  10  will  cause  a  progressive  shift  of  the  current  from  the  radial  gap  of  length 
L  to  the  region  near  the  free  surface.  Figures  10a  through  10c  present  contour  plots  for 
M  =  10  and  h0  =  -20,  -60,  and  -  100. 

The  nondimensional  fluid  velocity  profiles  for  M  =  2  and  M  =  10  for  negative  h0 
values  show  in  general  positive  flow,  in  the  positive  x-direction,  in  the  region  of  the 
channel  near  the  moving  wall  because  of  irreversible,  viscous  friction  effects  between  the 
wall  and  the  fluid.  Away  from  the  close  proximity  of  the  moving  wall  the  fluid  velocity 
profiles  are  negative,  in  the  minus  x-direction,  and  approach  zero  in  the  vicinity  of  the 

bottom  stationary  wall.  The  fluid  in  the  region  between  10<z<  13  and  0  <y  <  1  in 

— — ♦ 

general  moves  in  the  negative  x-direction  because  of  j  x  Ba  body  forces  caused  by  the 
currents  leaving  the  side  wall  to  the  tip  of  the  moving  wall  interacting  with  the  axial 
magnetic  field  B0,  and/or  the  viscous  interaction  with  the  fluid  near  the  center  region  of  the 
channel  moving  in  the  minus  x-direction.  The  nondimensional  fluid  velocity  moves  much 
faster  in  the  negative  x-direction  when  M  is  changed  from  2  to  10  because  of  the  increased 
j  x  B0  body  forces. 

The  nondimensional  induced  magnetic  field  contours  for  M  =  2  and  M  =  10  are  very 
similar  in  magnitude  and  position  for  the  same  negative  ha  except  in  the  side  region 
1 0  <  z  <  1 3  and  0  <  y  <  1 .  As  M  increases  from  M  =  2  to  1 0  for  the  same  negative  h0 
contour  lines  of  h(y,z)  leave  the  right  side  wall  and  approach  the  tip  of  the  top  moving  wall. 
In  each  figure  of  h(y,z)  for  M  =  2  there  is  one  contour  line  near  the  free  surface  that  leaves 
the  stationary  wall  and  moves  to  the  tip  of  the  moving  wall. 

Comparing  Figs.  5a  through  10c,  it  is  evident  that  the  results  for  M  =  2,  10  and 
h0  =  ±20,  ±40,  ±60,  ±80,  and  ±100  are  not  symmetrical  about  h0  =  0  for  the  same  M 
value.  For  h0  -  0,  jy  <0  for  0  <  z  <  9.5  and  jy  >0  for  9.5  £  z  <  13.  As  h0  increases 
from  0  to  100,  we  are  adding  a  positive  jy  everywhere.  Thus,  the  original  and  added 
current  density  cancel  for  0<z<9.5  and  reinforce  for  9.5  <  z  <  1 3.  As  h0  becomes 
negative  from  0  to  -  100,  we  are  adding  a  negative  jy  everywhere.  Thus  the  original  and 
added  current  density  reinforce  for  0  <  z  <  9.5  and  cancel  for  9.5  <  z  <  1 3.  As  a  result 
the  value  of\jy  \  for  0  <  z  <  9.5  is  larger  for  h0  =  -  20,  than  for  h0  =  +  20,  and  the  same 
M.  Similarily,  the  current  densities  near  the  free  surface  are  smaller  for  h,,  -  20  than  for 
h0  =  +  20  and  the  same  M. 

For  both  h0  =  ±20,  there  is  a  shift  of  current  density  from  the  radial  gap  to  the  free 
surface  region  as  M  increases  from  2  to  10.  However,  this  shift  is  not  as  strong  for  a  given 
change  in  M  for  h0  =  -20,  as  it  is  for  h0  -  +  20,  because  the  potential  difference  must 
overcome  the  tendency  for  current  to  flow  from  the  stationary'  to  the  moving  wall  near  the 
free  surface. 

These  same  general  principles  apply  for  the  higher  h0  values  at  the  same  M  value. 
Table  2  shows  for  comparison  the  values  calculated  for  simple  one-dimensional, 
nondimensional  Couette  flow  in  an  axial  magnetic  field  with  no  nondimensional  pressure 
gradient. 
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*  NONDIMENSIONAL  HEIGHT  y  NONDIMENSIONAL  HEIGHT  y  NONDIMENSIONAL  HEIGHT 


z  NONDIMENSION  At  WIDTH 


Fig.  10a.  M  =  10,  A  0  =  -  20 


Fig.  10b.  M  =  10,  h0  =  -  60 


Fig.  10c.  M  =  ’0 ,  h0  =  -  100 

Fig.  10.  Nondimensional  velocity  contours  u(y,z)  and  nondimensional 
current  density  stream  function  A(y,z)  for  Hartmann  number 
M  =  10  and  h0  =  -20,  -60,  and  -100. 
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Table  2.  One-dimensional  Couette  flow  values  for  axial  magnetic  field 


mm 

u(-75) 

u(.  50) 

w(.25) 

umin 

1 

-20 

+0.37 

0.0 

-0.13 

-0.13 

0.25 

-40 

-0.01 

-0.50 

-0.51 

-0.57 

0.38 

M  =  2 

-60 

-0.39 

-1.00 

-0.89 

-1.04 

0.42 

-80 

-0.77 

-1.50 

-1.27 

-1.53 

0.44 

-100 

-1.15 

-2.00 

-1.65 

-2.03 

0.45 

-20 

-1.15 

-2.00 

-1.65 

-2.03 

0.45 

-40 

-3.05 

-4.50 

-3.65 

-4.51 

0.48 

M=  10 

-60 

-4.95 

-7.00 

-5.55 

-7.01 

0.48 

-80 

-6.85 

-9.50 

-7.35 

-9.51 

0.49 

-100 

-8.75 

-12.00 

-9.25 

-12.01 

0.49 
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For  reference.  Appendix  C  presents  much  more  detailed  contour  plots  for  Figs.  5 
through  10.  The  figures  are  not  presented  as  composite  plots  of  the  contours  for 
nondimensional  velocity  and  nondimensional  current  density  stream  function,  as  they  are 
here,  but  as  single  plots  for  each  case. 

COMPARISON  OF  NUMERICAL  PARAMETERS  FOR 
ONE-DIMENSIONAL  AND  TWO-DIMENSIONAL  FLOW 

Tables  3  to  7  are  presented  for  various  numerical  parameters  for  both  one¬ 
dimensional  (see  Appendix  B)  and  two-dimensional  Couette  flow  in  an  external 
homogeneous  axial  magnetic  field  Ba  with  no  pressure  gradient  along  the  channel 

dp/dx  =  0. 

For  the  two-dimensional  case  the  nondimensional  length  of  the  perfectly  conducting, 
moving  wall  from  the  center  of  the  wall  to  the  end  is  b  =  10,  and  the  nondimensional  free 
surface  length  is  a  =  3,  (see  Fig.  2).  For  the  one-dimensional  case  the  perfectly 
conducting  moving  wall  nondimensional  length  from  the  center  of  the  wall  to  the  end  is 
b=  10.  We  have  compared  the  two  cases,  both  having  the  moving  wall  of  nondimensional 
length  b  =  10,  to  determine  the  end  effects  in  the  two-dimensional  case.  The  one¬ 
dimensional  case  has  no  end  effects. 

Each  of  Tables  3  to  7  is  for  a  specific  Hartmann  number  M  (where  M  can  run  through 
the  values  2,  4,  6,  8  and  10)  and  extends  through  a  range  of  h0  values  where  h0  =  -100, 

-  80,  -  60,  -  40,  -  20,  0,  20,  40,  60,  80,  100.  The  quantity  h0  is  a  measure  of  the  load 
current  (see  Eq.  16.2).  The  symbol  in  each  table  for  <p0  represents  the  two-dimensional, 
nondimensional  potential  at  the  perfectly  conducting  moving  wall  <p0  =  <p(y  =  l,z)  for 
-b<y<b  .  The  value  is  the  corresponding  one-dimensional,  nondimensional 
potential.  The  symbol  Pj  is  the  nondimensional  Joulean  dissipation  per  unit 
nondimensional  length  (Ax)  in  the  channel  for  the  two-dimensional  case.  Pjid  is  the 
corresponding  nondimensional  Joulean  dissipation  in  the  one-dimensional  case.  The 
symbol  P^  is  the  nondimensional  viscous  dissipation  in  the  two-dimensional  case  per  unit 
nondimensional  length,  and  P^xo  is  the  corresponding  one-dimensional  case. 

Comparison  of  the  numerical  data  for  the  one-dimensional  and  two-dimensional  cases 
in  Tables  3  to  7  leads  to  the  following  conclusions.  For  every  instance,  the  relationship 
I  <f>0 1  <  I  (p0]D  I  holds  for  the  nondimensional  potentials  in  the  two  cases.  The  one¬ 
dimensional  solution  assumes  that  all  current  flows  across  the  magnetic  field  lines  and  thus 
is  a  current  path  with  a  large  resistance.  In  the  two-dimensional  case,  a  fraction  of  the 
current  flows  along  the  magnetic  field  lines  near  the  free  surface  and  enters  the  moving  wall 
near  the  edge  at  z  =  b.  This  path  has  a  lower  resistance  so  that  the  overall  potential 
difference  is  smaller  than  the  one-dimensional  case.  Also,  the  fraction  of  the  total  current 
near  the  free  surface  increases  as  M  increases  so  I  <p0lD~<Po  I  increases  as  M  increases. 

Nondimensional  Joulean  power  losses  for  no  load  current  h0  =  0  are  Pjid  =  0  for 
the  one-dimensional  case  for  M  =  2,  4,  6,  8  and  10.  For  the  same  values  of  M,  Pj  extends 
from  0.28  through  -2.74  for  the  two-dimensional  case.  In  the  two-dimensional  case,  the 
weak  nondimensional  Joulean  power  losses  are  from  the  weak  eddy  current  when  there  is 
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Table  3.  Numerical  data  for  Hartmann  number  A/= 2 


mm 

5.962 

6.167 

1,802.1 

2,000 

650.9 

686.7 

mm 

4.674 

4.833 

1,153.5 

1,280 

424.1 

446.7 

-60 

3.385 

3.500 

o49.1 

720 

247.7 

260.0 

-40 

2.097 

2.167 

288.7 

320 

121.7 

126.7 

-20 

0.808 

0.833 

72.44 

80 

46.01 

46.67 

2 

0 

-0.481 

-0.500 

0.2793 

0 

20.71 

20.0 

2 

20 

1.769 

-1.833 

72.23 

80 

45.77 

46.67 

2 

40 

-3.058 

-3.167 

288.3 

320 

121.2 

126.7 

2 

60 

-4.346 

-4.500 

648.4 

720 

247.0 

260.0 

2 

80 

-5.635 

-5.833 

1,152.7 

1,280 

423.2 

446.7 

2 

100 

-6.924 

-7.167 

1,801.1 

2,000 

649.7 

686.7 

Table  4.  Numerical  data  for  Hartmann  number  M  =  4 


Hi 

5.072 

5.333 

1,830.2 

2,000 

2,488.6 

2,686.7 

Hi 

3.962 

4.167 

1,170.8 

1,280 

1,601.3 

1,726.7 

-60 

2.852 

3.000 

658.2 

720 

910.9 

980.0 

-40 

1.742 

1.833 

292.3 

320 

417.3 

446.7 

-20 

0.632 

0.667 

73.11 

80 

120.7 

126.7 

0 

-0.478 

-0.500 

0.719 

0 

20.92 

20.0 

20 

-1.587 

-1.667 

75.07 

80 

118.0 

126.7 

40 

-2.697 

-2.833 

296.2 

320 

412.1 

446.7 

60 

-3.807 

-4.000 

664.0 

720 

903.0 

980.0 

80 

-4.917 

-5.167 

1,178.6 

1,280 

1,590.8 

1,726.7 

100 

-6.027 

-6.333 

1,840.0 

2,000 

2,475.5 

2,686.7 

Table  5.  Numerical  data  for  Hartmann  number,  M  =  6 


M 

h0 

Pj 

P.JW 

Pu 

Pul  D 

6 

-100 

5.743 

6.167 

1,945.4 

2,000 

5,394.2 

6,020 

6 

-80 

4.500 

4.833 

1,242.2 

1,280 

3,463.6 

3,860 

6 

-60 

3.257 

3.500 

696.1 

720 

1,961.0 

2.180 

6 

-40 

2.014 

2.167 

307.3 

320 

886.4 

980 

6 

-20 

0.771 

0.833 

75.68 

80 

239.8 

260 

6 

0 

-0.472 

-0.500 

1 .257 

0 

21.164 

20 

6 

20 

-1.715 

-1.833 

84.04 

80 

230.6 

260 

6 

40 

-2.958 

-3.167 

324.0 

320 

868.0 

980 

6 

60 

-4.201 

-4.500 

721.2 

720 

1,933.4 

2,180 

6 

80 

-5.444 

-5.833 

1,275.6 

1,280 

3,426.8 

3,860 

6 

100 

-6.687 

-7.167 

1,987.2 

2,000 

5,348.2 

6,020 

Table  6.  Numerical  data  for  Hartmann  number  M  =  8 


M 

h0 

0(9 

_ 4 

)o\D 

P.1 

Pj\D 

.  Pn 

P  ulD 

8 

-100 

6.735 

7.417 

2239.9 

2000 

9162.8 

10,686.7 

8 

-  80 

5.295 

5.833 

1426.1 

1280 

5880.4 

6846.7 

8 

-60 

3.855 

4.250 

795.5 

720 

3325.0 

3860.0 

8 

-40 

2.415 

2.667 

347.9 

320 

1496.8 

1726.7 

8 

-20 

0.975 

1.083 

83.38 

80 

395.5 

446.7 

8 

0 

-0.465 

-0.500 

1.940 

0 

21.397 

20.0 

8 

20 

-1.905 

-2.083 

103.58 

80 

374.3 

446.7 

8 

40 

-3.345 

-3.667 

388.3 

320 

1454.3 

1726.7 

8 

60 

-4.785 

-5.250 

856.1 

720 

3261.4 

3860.0 

8 

80 

-6.225 

-6.833 

1506.9 

1280 

5795.5 

6846.7 

8 

-7.665 

-8.417 

2340.9 

2000 

9056.7 

10,686.7 

9 
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Table  7.  Numerical  data  for  Hartmann  number  M  =  10 


I 


M 

ho 

§0 

4>a1D 

P.! 

P.J\D 

Pu 

Pino 

10 

-100 

7.788 

8.833 

2,816.6 

2,000 

13,552.4 

16,686.7 

10 

-80 

6.139 

6.967 

1,788.3 

1,280 

8,697.1 

10,686.7 

10 

-60 

4.490 

5.100 

992.8 

720 

4,916.3 

6,020.0 

10 

-40 

2.842 

3.233 

430.1 

320 

2,210.2 

2,686.7 

10 

-20 

1.193 

1.367 

100.1 

80 

578.6 

686.7 

10 

0 

-0.456 

-0.500 

2.736 

0 

2,1.650 

20.0 

10 

20 

-2.104 

-2.367 

138.2 

80 

539.3 

686.7 

10 

40 

-3.753 

-4.233 

506.4 

320 

2,131.5 

2,686.7 

10 

60 

-5.401 

-6.100 

1,107.3 

720 

4,798.2 

6,020.0 

10 

80 

-7.050 

-7.967 

1,941.0 

1,280 

8,539.6 

10,686.7 

10 

100 

-8.699 

-9.833 

3,007.4 

2,000 

13,355.6 

16,686.7 
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no  load  current.  For  a  load  current,  h0  *  0  when  M  =  2,  4  and  6,  the  following 
relationship  is  found  to  be  true:  Pj  <Pjw-  The  one-dimensional  solution  assumes 
0  <  Iz  I  <  b,  but  in  the  two-dimensional  case  0  <  Iz  I  <  a  +  b,  which  allows  currents  to 
fringe  into  the  end  regions  giving  a  smaller  average  y-component  current  density  jy  and 
smaller  Pj.  For  hQ*  0  and  M  =  8  and  10,  the  relationship  Pjid<Pj  is  found  to  hold.  In 
the  two-dimensional  case  a  progressively  larger  fraction  of  current  is  concentrated  near  the 
free  surface  and  rotor  edge.  As  \j  I  becomes  locally  large  jy+jz  becomes  very  large 
and  this  extra  contribution  is  much  larger  than  the  reduction  elsewhere  in  the  channel. 

For  the  nondimensional  viscous  dissipation  when  the  load  current  is  not  zero  h0± 0, 
the  relationship  P^  <  P/xlD  is  found  to  hold.  The  one-dimensional  velocity  consists  of  the 
linear  velocity  y  plus  a  parabola  driven  by  the  body  force  jyB0.  In  P^d  (Appendix  B, 

Eq.  B-l  1)  the  2b  is  due  to  the  linear  velocity  u  =  y  and  M2  is  due  to  the  parabola. 

In  the  two-dimensional  case  as  M  increases,  an  increasing  fraction  of  the  current  flows 
along  the  field  lines  near  the  free  surface.  This  current  produces  no  velocity  except  near  the 
moving  wall  edge  where  it  must  cross  field  lines  to  enter  the  wall.  Thus  in  the  two- 
dimensional  case,  the  parabola  near  z  =  0  is  much  smaller  than  in  the  one-dimensional 
case,  and  P ^  is  concentrated  near  the  edge  of  the  moving  wall. 

The  one-dimensional  variables  calculated  and  presented  in  this  section  were  derived  in 
Appendix  B,  but  are  presented  here  for  continuity  of  presentation.  The  nondimensional  y- 
component  of  current  density  is 

- 

,yW  ~  Mb  ’  (26) 


and  the  nondimensional  Joulean  power  loss  per  unit  nondimensional  length  is  given  by 
p _ 2 h} 


The  nondimensional  velocity  is 

and  the  nondimensional  viscous  power  loss  per  unit  nondimensional  length  is  given  by 
M2h2 

P,«n  =  2b  +  ^-J^ 


The  nondimensional  potential  in  the  channel  has  the  following  form 

y2  IMh0\jy2  y3  I  h0  \ 

2  \  2b  2  ~  3  ~\Mb\y  ’ 

and  the  nondimensional  potential  at  the  moving  wall  is 

^id  =  -“(i^)(t7  +  7t)  -  where  (<Pid  (>’  =  0  =  <Po1d)- 


2  \bj\M  12/ . 

(Note:  all  the  variables  in  these  equations  were  defined  earlier  in  the  report.) 


(30) 

(31) 
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CONCLUSIONS 


Numerical  data  were  presented  for  the  nondimensional  velocity  profiles  u(y,z)  and 
nondimensional  induced  magnetic  field  h  (y,z)  for  an  incompressible  conducting  fluid  in  a 
two-dimensional  magnetohydrodynamic  channel  in  an  axial  homogeneous,  external 
magnetic  field  (Couette  flow).  The  top  wall  of  the  channel  was  a  moving,  perfectly 
conducting  wall  with  insulating  free  surfaces  in  both  top  comers.  The  bottom  stationary 
wall  and  the  right  and  left  stationary  sidewalls  were  also  perfect  conductors.  A  current  load 
was  simulated  between  the  top  and  bottom  conducting  walls  by  applying  an  external 
potential  to  the  top  wall.  Three  cases  for  Hartmann  numbers  M  =  2  and  10  were  studied 
for  a  load  current  moving  from  bottom  to  top  wall  (h0  =  20,  60,  100),  top  to  bottom  wall 
(h0  =  -20,  -60,  - 100)  and  no  load  current  ( h0  =  0).  The  data  were  calculated  for  a 
rectangular  channel  where  the  radial  gap  was  chosen  as  length  L ,  each  axial  free  surface 
gap  as  3 L  ,  and  the  length  of  the  moving  wall  as  20L.  In  this  work  the  nondimensional 

pressure  gradient  ^/dx  was  set  equal  to  zero.  Thus,  the  fluid  motion  was  due  only  to  the 
drag  of  the  moving  wall  and  the  Lorentz  body  force. 

These  three  sets  of  data  supported  the  concept  that  the  radial  gap  and  the  tree  surface 
region  represent  electrical  resistances  in  parallel  between  the  stationary,  perfectly 
conducting  walls  and  the  moving,  perfectly  conducting  wall.  The  data  also  showed  that  the 
resistance  of  the  radial  gap  increases  as  M  2  while  that  of  the  free  surface  region  increases 
by  approximately  Thus  as  M  increases,  the  division  of  current  shifts  to  the  free 
surface  region  and  the  current  density  in  the  radial  gap  decreases  as  M~l. 

Numerical  data  for  one-dimensional  Couette  flow  with  no  pressure  gradient  in  an 
axial  magnetic  field  having  no  end  effects  were  presented  and  compared  with  numerical 
data  for  two-dimensional  Couette  flow. 

On  the  basis  of  the  theory  and  the  range  of  data  computed,  a  number  of  conclusions 
can  be  drawn: 

1.  For  every  case  studied  the  absolute  value  of  the  nondimensional  potential  at  the 
perfectly  conducting  moving  wall  of  the  two-dimensional  case  I <p0\  is  less  than  the 
corresponding  nondimensional  potential  \<p0JD I  in  the  one-dimensional  case.  The 
fraction  of  the  total  current  near  the  free  surface  in  the  two-dimensional  case 
increases  as  M  increases,  so  \<p0lD  ~  <Po\  increases  as  M  increases. 

2 .  For  the  range  of  data  computed,  with  no  load  current  across  the  channel,  h0  =  0, 
the  nondimensional  Joulean  power  losses  per  nondimensional  length  Pjjd  in  the 
one-dimensional  case  are  zero.  The  corresponding  Joulean  losses  per 
nondimensional  length  Pj  for  the  two-dimensional  case  extend  from  0.28  at 

M  =  2  to  2.74  for  M  =  10.  The  power  losses  in  the  two-dimensional  case  are 
caused  by  the  weak  eddy  currents.  For  load  currents  across  the  channel,  the 
power  relationship  Pj  <Pj\d  holds  when  ho*0\  for  M  =  2,  4  and  6  and  Pjjd 
<Pj  holds  for  M  =  8  and  10. 

3 .  It  was  found  that  the  nondimensional  viscous  dissipation  per  nondimensional 
length  Pjj  is  always  less  than  the  corresponding  one-dimensional  case  for  the 
range  of  data  computed. 
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These  results  show  that  it  is  necessary  to  model  the  magnetohydrodynamic  flow  in  a 
current  collector  by  a  two-dimensional  model,  as  opposed  to  a  one-dimensional  model,  to 
obtain  a  fundamental  understanding  of  the  flow  profile,  current  density  distribution  in  the 
liquid  metal,  and  viscous  and  Joulean  power  losses. 
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APPENDIX  A 

MATHEMATICAL  SOLUTION  FOR  RECTANGULAR  CHANNEL 
WITH  FOUR  PERFECTLY  CONDUCTING  CONNECTED  WALLS, 

ONE  A  MOVING  WALL 

The  problem  discussed  herein  consists  of  a  rectangular  channel  in  a  homogeneous, 
external,  transverse  magnetic  field  B0  with  four  thin  perfectly  conducting  walls  (see 
Fig.  A-l).  The  perfectly  conducting  wall  at  the  top  of  the  figure  is  assumed  to  move  with 
a  velocity  u  =  u(z  =  za).  It  should  be  carefully  noted  that  the  coordinate  system  is  in  a 
different  position  in  regard  to  the  rectangular  channel  than  is  the  coordinate  system  in  the 
text  of  this  report. 

The  mathematical  problem  developed  in  this  Appendix  formed  a  basis  for  some  of  the 
theoretical  development  of  the  more  complex  model  presented  in  the  text  of  the  report.  The 
major  fault  of  the  mathematical  model  presented  herein  is  that  the  four  perfectly  conducting 
walls  are  connected  so  that  there  are  no  insulating,  free  surface  gaps  between  the  perfectly 
conducting  moving  wall  and  the  perfectly  conducting  left  and  right  side  stationary  walls  as 
are  in  the  rectangular  channel  mathematical  model  presented  in  the  text  of  the  report.  Thus 
in  the  model  in  this  Appendix,  it  is  impossible  to  apply  an  external  potential  <p0  to  the 
moving  walls  of  the  channel.  This  results  in  a  severe  limitation  to  the  mathematical  model 
since  current  collectors  carry  a  heavy  load  current. 

We  were  able  to  obtain  series  solutions  for  the  flow  velocity,  induced  magnetic  field, 
total  current,  and  voltages  in  the  channel.  These  solutions  are  presented  in  the  Appendix  to 
provide  a  continuity  and  theoretical  background  to  the  development  of  the  theoretical 
mathematical  model  presented  in  the  text  of  the  report.  Numerical  computer  results  for  the 
problem  herein  were  not  obtained. 

The  mathematical  problem  is  made  more  tractable  by  transforming  the  two 
dimensional,  steady,  incompressible,  viscous  magnetohydrodynamic  flow  partial 


differential  equations  into  nondimensional  form  using  the  notation  of  the  Hughes  and 
Young.10  They  worked  the  similar  rectangular  channel  problem  where  all  four  connected, 
perfectly  conducting  walls  are  stationary. 

The  nondimensional  variables  are  defined  as 

y0 

y2o  (sp/sx) 

(A- 1.1) 

V; 

• 

II 

If-  ^ 

yl(3pL 

(A-l. 2) 

z*  -  2 

(A-l. 3) 
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Fig.  A-l.  Channel  with  four  perfectly  conducting  walls. 
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where  the  dimensional  variables  (without  superscript  *)  have  the  following  meaning: 

y0  =  one  half  width  of  channel,  (A- 1 .4) 

z0  =  one  half  height  of  channel,  (A- 1.5) 

k  =  z°/y0  =  aspect  ratio,  (A- 1.6) 

jU/  =  viscosity  of  fluid  in  channel,  (A- 1 .7) 

dp/dx  =  constant  pressure  gradient  in  x-direction,  A- 1 .8) 

<T  =  conductivity  of  fluid,  (A- 1.9) 

B0  =  magnitude  of  the  homogeneous  external  magnetic  field,  A- 1 . 1 0) 

u  =  dimensional  fluid  velocity,  (A- 1.11) 

Hx  =  dimensional  induced  magnetic  field.  (A- 1.12) 

The  nondimensional  magnetohydrodynamic  equations  for  it*(y*,z*)  and  h*(y*,z *)  are 


H*x  d2H*x  tMdu' 
—  + - -  +  M  — 

>*2  dz*2  dy* 


(A-2.1) 


a  V  a v  3h*x  , 

av*2  dz*2  av* 


(A-2.2) 


where  the  appropriate  boundary  values  must  be  satisfied.  The  well-known  assumed 
nondimensional  solutions  have  the  following  form 


u*  =  it*  +  ^  Un{z*)  cos\ny* 
k—yo°  n  =  o 


(A-3.1) 


Hl  =  H*x  +  X  H*n(z*)sinX*y*  , 

£  — >  oo  n  =  0 

where 

ul{z*)  =  C\  sinhP*z*  +  C2 coshP*2z*  +  C->,sinhP\z*  +  C4  coshPlz*  , 

pi  pi  =  V  a;2  ±jC  m 


P\,P'2  =r*n±js*n 

,  *  _  (In  +  1)tt 
An  ~  2 

*  -  1  IcoshMy* 

M 2  '  coshM 


Hi 

k— 


sinhMy* 
M  coshM 


(A-3.2) 

(A-3.3) 

(A-3.4) 

(A-3.5) 

(A-3.6) 

(A-3.7) 

(A-3.8) 
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The  particular  solutions,  u  and  Hx  ,  are  the  solution  to  the  Row  problem 

k  -*  ~  k-*<x> 

between  two  perfect  conducting  infinite  plates.  Equations  A-3.1  and  A-3.2,  transform 
Eqs.  A-2.1  and  A-2.2  into  ordinary  differential  equations. 

The  solution  for  un*  is  assumed  unsymmetrical  with  respect  to  z*  and  was  thus 
expressed  as  a  combination  of  even  and  odd  functions 


u*  =  C\  cosh (r*  +  js*n)z *  +  C2  cosh (r*  - /y*)z *  + 

C3  sinh (r*  +  js*)z *  +  C4  sinh (r*  -  js* )z *  , 

or  simply  as 

«*(z * )  =  Ansinh  (r*  z * )  z  * )  +  Bn  cosh  (r*  z * )  co^s*n  z * ) 

Cnsinh (r* z * )  cos[sn  z * |  +  Dn cosh (r*  z  * )  z *  j 


(A-4.1) 


(A-4.2) 


where 


An=j(C i  -  C2)  , 
Bn  =  (C\  +  C2)  , 
Cn  =  (C3  +  C4)  , 
h>n  ~  j  (C3  —  C/\)  . 


(A-4.3) 


The  solution  for  Hn*(z  *)  is  also  assumed  unsymmetrical  with  respect  to  z  *  and  was 
expressed  as 


MXn 


where 


or  more  simply  as 


Xn  )ci  coshP*z*  +  [p^2  -  Xn  jc2coshP2z* 

Xn  )C\COShP*Z*  +[p22  -  Xn  )c2COShP2Z* 

(A-5.1) 

'*  -  j  C2  coshP^  z*+jC3  sinhP*  z *  -j  C4  sinhP 2  z *] 

(A-5.2) 

and  P*22  -  Xn  =  -j[xnM\  , 

(A-5.3) 

)  sin  z *)  —  An  cosh  (r^  z*)  cos  (s„  z*) 

+  Cncosh  (r„z*)sm  (j*  z*)  -  D„sinh(r^z*)  cos  (s^z*)  (A-5.4) 

The  following  boundary  conditions  must  be  satisfied  for  this  specific  problem 

00 

1  =  u*  +  wn(z*=  k)  cosXny*  , 


i  — ><»  «=o 


(A-6.1) 


1 
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(A-6.2) 


0  =  u*  +  X  USZ*~  -  k)  cosXny*  , 


n=0 


Q  dH*(y\z*=lc)_  x*  dHn(z*=  k) 
dz *  „=o  "  dz* 


0  _  dH*(y*,z*=  -k)  _  y  dHn\z'=-k) 
dz*  h  n  dz * 


cosXny *  , 

* 

cosXny * 


(A-6.3) 

(A-6.4) 


Equations  A-6.1  and  A-6.2  show  the  no-slip  boundary  conditions  on  the  conducting 
fluid  that  must  be  satisfied  along  the  moving  wail  at  z  *  =  k  and  the  stationary  wall  at 
z*  =-k.  Equations  A-6.3  and  A-6.4  show  that  the  y-component  of  the  current  density 
must  be  zero  along  these  perfectly  conducting  walls  at  z  *  =  k  and  -k  ,  since  it  is  well 
known  that  for  a  perfect  conductor  the  tangential  component  of  the  current  density  must  be 
zero  or  infinite  currents  will  flow  in  the  conductor. 

The  form  of  the  mathematical  solution  for  u*(y*,z*)  and  jz*(y*,z*)  guarantees  that 

u*(y*=  ±1,  z*)  =  0,  and 


/z(y*=±l,z*)  =  0 


(A-7.1) 

(A-7.2) 


on  both  the  left  and  right  wall.  Substitution  of  the  expressions  for 
u*  ,  Un,  H*x  ,  and  H*n,  derived  previously  in  this  report,  into  the  boundary  conditions  of 

k— >oo 


k — > 


Eqs.  A-6.1  to  A-6.4  determines  the  coefficients,  which  are 


An  -  En 

* 

rn  cosh 

r> 

1  sin 

?nk 

+  sn  sink 

r» 

cos 

Bn  =  En 

rnsinh\ 

cos 

S*k 

* 

-  sn  cosh 

K* 

i  sin 

s'A 

Cn  =  En 

r*n  cosh 

|  cos 

(*»* 

)  -  sn  sinh 

K* 

sin 

Dn  -  F  n 

r*nsinh\ 

sin  ( 

Snk ) 

+  sn  cosh 

cos 

'FA 

(A-8.1) 

(A-8.2) 

(A-8.3) 

(A-8.4) 


where 


En  = 


2(- 

\r 

„  .21 
2  +  M  +  Xn 

K 

M1  +  A*2) 

* 

rn 

sinh 

-  s*n  sin 

Fn  = 


2(-i  r 


Xn[r*n  sinh  (2 /•«£)  +  s„  sin  (2j^£)] 
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APPENDIX  B 

ONE-DIMENSIONAL  COUETTE  FLOW 
WITH  AN  AXIAL  MAGNETIC  FIELD, 

APPLIED  EXTERNAL  POTENTIAL,  AND  NO  PRESSURE  GRADIENT 

The  problem  discussed  herein  is  not  commonly  reported  in  the  literature,  and  is 
presented  here  for  continuity  of  presentation.  One-dimensional  Couette  flow  with  an 
externally  applied  potential  in  a  constant  homogeneous  axial  magnetic  field  and  no  pressure 
gradient  is  shown  in  Fig.  B-l.  The  nondimensional  externally  applied  potential  of  the 
moving,  perfectly  conducting  infinite  plate  is  < p0,  and  the  nondimensional  potential  of  the 
stationary,  perfectly  conducting  infinite  plate  is_zero  (0  =  0).  The  axial,  constant, 
homogeneous  externally  applied  magnetic  field  B  is  in  the  z-direction.  The  conductivity 
and  viscosity  of  the  liquid  metal  are  o  and  fJf,  respectively.  The  nondimensional  velocity  of 
the  liquid  metal  u(y)  is  in  the  ^-direction.  There  are  no  plate  end  effects  on  the  fluid  flow 

0 

since  the  plates  are  of  infinite  length  ( i.e.  —  =  0  ,  and  j2  =  0). 

dz 


y 


Fig.  B-l.  One-dimensional  Couette  flow  with  an  axial  magnetic  field 


The  numerical  data  calculated  for  the  variables  of  interest  for  one-dimensional  Couette 
flow  presented  here  will  be  compared  with  the  corresponding  numerical  data  for  the  two- 
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dimensional  Couette  flow  presented  in  the  text  of  the  report.  End  effects  are  important  in 
the  two-dimensional  Couette  flow  because  of  the  gaps  between  the  moving,  perfectly 
conducting  finite  top  plate  and  the  perfectly  conducting  side  walls.  The  two  gaps  are 
insulating,  stress-free,  free  surfaces.  In  the  one-dimensional  Couette  flow,  end  effects  are 
not  incorporated  in  the  derivation.  Therefore,  it  would  be  expected  that  the  two  sets  of  data 
should  be  similar  only  at  low  Hartmann  numbers  M. 

Referring  to  the  two-dimensional  magnetohydrodynamic  equations  in  the  text  of  the 
report,  the  governing  system  of  nondimensional  ordinary  differential  equations  for  one¬ 
dimensional  Couette  flow  with  a  constant,  homogeneous  axial  magnetic  field  is 


M~ 2 


d2u(y) 
dy 2 


+jy{y)  =  Q> 


jy{y)  = 


-d<p(y) 

dy 


-  u{y) . 


dJy(y)  _  0 

dy 


(B-l.l) 

(B-1.2) 

(B-1.3) 


The  boundary  conditions  for  the  nondimensional  velocity  u  and  nondimensional 
potential  <p  for  the  moving,  perfectly  conducting  top  wall  are 


uiy-  1)  =  1  , 


<P(y=  1)  =  , 


(B-2.1) 


and  the  stationary,  perfectly  conducting  bottom  wall  are 

U(y=0)  =  0,  0(y=O)  =  O.  (B-2.2) 


The  z-component  of  the  current  density  jz  is  equal  to  zero  (jz  =  0). 
From  Eqs.  B-1.2  and  B-1.3,  it  is  noted  that 


d2(p(y)  _  du(y) 
dy2  dy 


(B-2.3) 


The  mathematical  solutions  to  the  above  Couette  problem  for  the  liquid  metal 
nondimensional  velocity  u(y)  the  nondimensional  potential  in  the  fluid  <p(y ),  and  the  y- 
component  of  the  nondimensional  current  density  jy  are 


u(y)=y- 


1 + 12A/~2 


(B-4.1) 


0(y)  =  - 


2 


)  * 


(B-4.2) 


/  120o  +  6 1 

V 12  +  M2I 


(B-4.3) 


where  h0  is  defined  in  Eq.  B-7. 


40 


A  magnetic  stream  function  can  be  defined  by  h(z)  where  h(z)  is  a  linear  function  of  z 
(i.e.  constant  multiplied  by  z).  The  jy  and  jz  can  be  expressed  as 


_  _  6(20O+  l)  j  3/t(z)  h^_ 

y  M2  +  12  dz  Mb 


jz  =  0  =  M 


_i  Bh(z) 
dy 


(B-5.1) 

(B-5.2) 


Then  /i(z)  has  the  following  form 
6(2<^+  1 )  z 


h,= 


M  +  l2Af‘ 


-l 


(B-6) 


The  total  dimensionless  current  from  the  stationary  wall  to  the  moving  wall  is 

rb 


I 


=  l(Ax)j  jy(\ ,z)  dz  =  2{Ax)  M  ~xh0  . 


(B-7) 


Therefore  h0  can  be  expressed  as 
6  (2 (p0  +  l)  z 


h{z)  =  ~ 


M+  12M 


(B-8) 


where  b  is  the  length  of  a  finite  plate  in  the  z-direction  from  the  center.  From  Eq.  B  8  the 
nondimensional  potential  of  the  top  plate  <p0(y  =\)  =  (p0is 


2  1  lb 


(B-9) 


The  nondimensional  velocity  u(y )  in  terms  of  ha  is 

“(>’)  =  y  +  y  ( 1  -  y)  > 


(B-10) 


Using  Eq.  B-10,  the  total  viscous  dissipation  per  nondimensional  unit  length  in  the 
x-direction  for  plates  of  length  2b  is 


Pn  ID 


-J7 


^dydz-lb^l 


dy  ) 


6  b 


(B-ll) 


ID 


where  P„  w  = 

P/Uo 
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Also  using  Eq.  B-5,  the  total  Joulean  dissipation  per  nondimensional  unit  length  in  the  x- 
direction  for  plates  of  length  2b  is 


Pjw=M2 


jy  dy  dz 


(B-12) 


where  Pj\o  = 


and  the  superscript  (*)  represents  a  dimensional  quantity. 

To  understand  the  relationship  between  the  nondimensional  applied  external  potential 
at  the  moving  wall  (p0  and  h0  for  the  one-dimensional  Couette  solution,  we  must  remember 
that  a  constant  h0  represents  a  constant  total  current,  but  a  constant  <p(y  =  1)  =  <p0 
represents  a  potential  difference  that  increases  linearly  with  B0, 

<Po  =  U0B0L<p0  (B-13) 


hn  = 


1* 


2£/*V^(Ax*)  ’ 


where  I*  =  total  dimensional  current  from  stator  to  rotor  over  rotor  length  2bL.  The 
supercript  (*)  represents  a  dimensional  quantity.  The  dimensional  applied  potential  to  the 
moving  wall  <p0*  is 


<Po  =-'/2  U0B0L- 


2b  (Ax*)  a 


1  + 


M2 

12 


(B-14) 


Here  -  xh  U0  B0  L  is  the  potential  difference  needed  to  balance  the  average  of  u  x  B0  if  u 
is  simple  Couette  flow,  which  it  is  for/*  =0  in  the  one-dimensional  solution. 


The  resistance  of  the  stagnant  liquid  metal  for  plates  of  length  2bL  is 


2b  (Ax*)  o 


(B-15) 


When  M  1  ,  the  current  does  not  change  the  velocity  profile,  so 
(j)*o  +  '/2  U0B0L  =  -  I*  R*  . 


(B-16) 
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As  the  Hartmann  number  M  increases,  the  electric  current  produces  a  body 
force  jy  B0x  which  accelerates  the  flow  so  that  u  deviates  from  the  simple  Couette  flow 
(u(y)=y).  The  change  in  velocity  always  produces  an  induced  electric  field  - u*B0y , 
which  opposes  the  current.  If  jy>  0,  then  -u*B0y  is  negative;  if jy<0,  then  -u*Bay  is 
positive.  Therefore  the  moving  wall  and  the  stationary  wall  see  a  resistance  which  is  more 
than  that  of  stagnant  liquid  metal 

1+^-1.  (B- 17) 


Reff=R* 
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We  will  see  that  R*ft  actually  increases  roughly  as  M  because  of  end  currents  in 
thetwo-dimensional  Couette  problem  in  the  text  of  the  report.  The  expression  for  the 
nondimensional  velocity 


u(y)  =  y  + 


Mho 

2b 


y(i  - y ) 


(B-18) 


consists  of  the  linear  Couette  term  y,  plus  a  parabolic  term.  The  maximum  and  minimum 
points  for  u(y),  um  are  given  by 


Mh0 

4b 


(B- 1 9.1) 


.  1  b 

where  ym  =  ^  +  jj7~ 

2  A i  Hq 


(B-19.2) 


and  M>-~  for  ho<0,  M>^~  for  ho>0 
h0  h0 


(B-19.3) 


where  the  subscript  m  represents  a  maximum  or  minumum  point.  For  ho<0,M  must 
have  a  certain  minimum  value  before  a  negative  velocity  will  occur.  For  h0  >  0,  M  must 
have  a  minimum  before  umax  is  greater  than  1 .  Two  typical  plots  of  u(y)  versus  y  are 
shown  in  Fig.  B-2  for  h0>  0  and  h0  <0.  The  values  of  M  are  high  enough  in  both  cases 
for  a  maximum  umax  for  h0  >  0,  and  minimum  umin  for  h0  <  0. 


Typical  values  of  Couette  flow  in  an  axial  magnetic  field  with  no  pressure  gradient 
dh are  presented  in  Tables  B-l  and  B-2.  Figure  B-3  presents  two  non-dimensional 


velocity  profiles  w(y)  versus  nondimensional  distance  between  infinite  plates  y  for  one¬ 
dimensional  Couette  flow  for  h0  =  ±  10  and  Hartmann  number  M  =  1.  In  this  case,  the 
Hartmann  number  M  is  so  low  that  there  are  no  maximum  or  minimum  points  on  the 
profiles. 
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/min  y=  1 

NONDIMENSIONAL  DISTANCE  BETWEEN  PLATES 


Fig.  B-2.  Typical  nondimensional  velocity  profiles  for  h0>{)  and  h0  <  0. 
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Table  B-l.  One-dimensional  Couette  flow  velocity  for  b  =  10 


NONDIMENSIONAL  VELOCITY  [u(y)] 


APPENDIX  C 

COMPLEX  CONTOUR  PLOTS 


This  Appendix  presents  the  complex  contour  plot  of  the  nondimensional  velocity 
profile  for  a  specific  M  and  h0  at  the  top  of  each  page.  At  the  bottom  of  the  page  is  the 
corresponding  contour  plot  of  the  current  density  stream  function.  The  plots  are  sorted 
according  to  the  following  5  groups: 


> 

Figures  C-1A  -  C-2b 

M  =  2, 

ho  ~ 

0 

M  =  10, 

ho  = 

0 

> 

Figures  C-3a  -  C-7b 

M  =  2, 

he  = 

20 

M  =  2, 

hQ  - 

40 

M  =  2, 

ho  = 

60 

M  =2, 

h0  = 

80 

M  =  2, 

h0  = 

100 

> 

Figures  C-8a  -  C-  12b 

M  =  10, 

ho  = 

20 

M  =  10, 

ho  — 

40 

M  =  10, 

ho  ~ 

60 

M  =  10, 

h0  = 

80 

M  =  10, 

ho  = 

100 

> 

Figures  C-  13a  -  C-  17b 

M  =  2, 

h0  — 

-20 

M  =  2, 

ho  ~ 

-40 

M  =  2, 

ho  ~ 

-60 

M  =  2, 

h0  - 

-80 

M  =  2, 

h0  = 

-  100 

> 

Figures  C- 1 8a  -  C-22b 

M  =  10, 

ho  = 

-20 

M  =  10, 

ho  = 

-40 

Af  =  10, 

ho  — 

-60 

M  =  10, 

h0  = 

-80 

M  =  10, 

ho  = 

-  100. 

The  symbol  M  represents  the  Hartmann  number,  and  h0  is  a  measure  of  the  load  current 
across  the  rectangular  channel. 
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y  NONDIMENSIONAL  HEIGHT  y  NONDIMENSIONAL  HEIGHT 
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Fig.  C-2a.  Nondimensional  velocity  contours  u(j,z)  for 
Hartmann  number  M  =  10  and  h0  =  0. 
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Fig.  C-2b.  Nondimensional  current  density  stream  function  contours 
h(y,z)  x  103  for  Hartmann  number  M  =  10  and  h0  =  0. 
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Fig.  C-3a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  2  and  h0  =  20. 
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Fig.  C -31).  Nondimensional  current  density  stream  function  contours 
h(y,z )  for  Hartmann  number  M  =  2  and  h0  =  20. 
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y  NONDIMENSIONAL  HEIGHT  y^lONDIMENSIONAL  HEIGHT 
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Fig.  C-4a.  Nondimensional  velocity  contours  u(y,z)  f»r 
Hartmann  number  M  =  2  and  h0  =  40. 
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Fig.  C-4b.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  2  and  h0  =  40. 
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Fig.  C-5a.  Nondimensional  velocity  contours  u (>•,:)  for 
Hartmann  number  M  -  2  and  h0  -  60. 
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Fig.  C-5b.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  -  2  and  h0  =  60. 
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Fig.  C-6a.  Nondimensional  velocity  contours  «(y,z)for 
Hartmann  number  M  =  2  and  h0  =  80. 
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Fig.  C-6b.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  2  and  h0  -  80. 
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Fig.  C-7a.  Nondimcnsional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =2  and  h0  =  100. 
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Fig.  C-71).  Nondimcnsional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  -  2  and  h0  -  100. 
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Fig.  C-8a.  Nondimensionai  velocity  contours  u(j,z )  for 
Hartmann  number  M  =  10  and  h0  =  20. 
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Fig.  C-8b.  Nondimensionai  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  Af  =  10  and  h0  =  20. 
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Fig.  C-9a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  10  and  h0  =  40. 
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Fig.  C-9b.  Nondimensional  current  density  stream  function  contours 
li(y,z)  for  Hartmann  number  M  -  10  and  hn  -  40. 
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Fig.  C-lOa.  Nondimensional  velocity  contours  u(y,z) 

for  Hartmann  number  M  =  10  and  h0  =  60. 


I 

0.4  S  ! 

I  l 


0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0  11.0  12.0  13.0 

z  NONDIMENSIONAL  WIDTH 

Fig.  C-lOb.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  10  and  h0  =  60. 
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Fig.  C-lla.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  10  and  h0  =  80. 
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Fig.  C-llb.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  10  and  h0  =  80. 
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l-ig.  CM 2a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  -  10  and  h0  =  100. 
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l  i(-.  (’-12b.  Nondimensional  current  density  stream  function  contours 

IHy,z)  for  Hartmann  number  M  -  10  and  h0  -  100. 
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Fig.  C-13a.  Nondimensional  velocity  contours  «(y,z)  for 
Hartmann  number  M  =  2  and  h0  -  -  20 . 
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Fig.  (  -131).  Nondimensional  current  density  stream  function  contours 
h(y.z )  for  Hartmann  number  M  =2  and  h0  -  -20. 
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Fig.  C-14a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  2  and  h0  =-40. 
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Fig.  C-141).  Nondimensional  current  density  stream  function  contours 
li(y,z)  for  Hartmann  number  M  =  2  and  h0  =-40. 
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Fij>.  C-15a.  Nondinunsional  velocity  contours  u(y,z)  for 
Hartmann  number  M  -  2  and  h0  =-60. 
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l  ift.  C-151).  Nondinunsional  current  density  stream  function  contours 
/i(>',z)  for  Hartmann  number  M  =  2  and  h0  -  -60. 
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Fig.  C-16a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  2  and  h0  =  -80. 
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Fig.  C-16t>.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  2  and  li0  =-80. 
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Fig.  C-17a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  2  and  h0  =  -  100. 


z  NONDIMENSIONAL  WIDTH 

Fig.  C-17b.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  2  and  Ao  =  -100. 
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Fig.  C-18a.  Nondimensional  velocity  contours  u[y,z)  for 
Hartmann  number  M  =  10  and  h0  -  -20. 


Fig.  C-18b.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  10  and  li0  =  -20. 
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Fig.  C-19a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  10  and  h0  =  -40. 
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Fig.  C-191).  Nondimensional  current  density  stream  function  contours 
hty,:)  for  Hartmann  number  M  =  10  and  h0  =  -40. 
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Fig.  C-20a.  Nondimensional  velocity  contours  u{y,z)  for 
Hartmann  number  M  =  10  and  h0  =  -60. 
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Fig.  C-201).  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  Af  =  10  and  h0  -  -60. 


67 


v  NONDIMENSIONAL  HEIGHT  V  NONDIMENSIONAL  HEIGHT 


MOVING  WALL,  u  =  1 


FREE  SURFACE -H 


1.0  . 

1 

1 

i  :  i  '  0.0  1  ;  -  : 

T  .  «  0.0  :  • 

! 

-2.40 

-2.40 

0.8 

-4.80 

-4.80 

1 

"7.20 

-7.20 

-7.20 

0.6 

§ 

0.4 

V 

/ 

-7.20 

-7.20 

<v 

/ 

0.2 

t 

-4.80 

-4.80 

j 

r  -2.40 

-2.40 

0.0 

. -1 - v.  .  — -  —  -  .1  - 

1  .  1  i 

0.0 

0.0  1.0  2.0  3.0  4.0 

5.0  6.0  7.0  8.0 

9.0  10.0 

11.0  12.0 

z  NONDIMENSIONAL  WIDTH 

13.0 


Fig.  C-21a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  10  and  h0  =  -80. 
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F'ig.  C-21b.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  10  and  h0  =  -80. 
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Fig.  C-22a.  Nondimensional  velocity  contours  u(y,z)  for 
Hartmann  number  M  =  10  and  h0  =  -  100. 
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Fig.  C-22b.  Nondimensional  current  density  stream  function  contours 
h(y,z)  for  Hartmann  number  M  =  10  and  h0  =  -  100. 
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